#ifndef POLY_DEF_H
#define POLY_DEF_H
/*
 * include this file in 
 *  poly\balance.c
 *  poly\qr.c
 */
#include <stdlib.h>
#include "gsl\gsl_complex.h"
#include "gsl\gsl_errno.h"
#include "gsl\gsl_machine.h"

/* C-style matrix elements */
#define MAT(m,i,j,n) ((m)[(i)*(n) + (j)])

/* Fortran-style matrix elements */
#define FMAT(m,i,j,n) ((m)[((i)-1)*(n) + ((j)-1)])

#endif

